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ABSTRACT 


An experiment was conducted in the Florida Keys from 2 - 19 July 1990 to test 
the performance of the Naval Oceanic Vertical Acrosol Model (NOVAM) in a weak- 
convective regime. Meteorological data collected by aircraft and boat was used to gen- 
erate the surface and vertical profile information files required by NOVAM. Using this 
information, NOVAM predicts the aerosol extinction (km!) for a vertical cross-section 
of the atmosphere. Arrcraft-observed aerosol extinction profiles were also obtained. 
Comparisons between observed and NOVAM aerosol extinction profiles revealed major 
deviations above the cloud top. From the surface to the top of the cloud layer, NOVAM 
generally did an excellent job in predicting profile shape, with the magnitude of aerosol 
extinction tied to the extinction matched at the surface. In a few cases, observed ex- 
tinction increased more rapidly than NOVAM predicted extinction from the surface to 
the base of the cloud laver. This is attributed to ram scavenging associated with 
thunderstorm activity in the area. Comparison between different aerosol extinction 
profiles revealed much spatial and temporal variation that was verified by Lidar profiles 
of atmospheric structure. Thunderstorm activity, multiple cloud-lavers, and the spatial 
variation in the atmospheric structure have led to a hvpothesis that deep-convection was 
responsible for the major differences between observed and predicted aerosol extinction 
profiles. If this is the case, a simple modification to the weak-convective model may be 


made to apply to a deep-convection model. 
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I. INTRODUCTION 


Successful employment of modern weapon systems require the continuous ability to 
use wavelengths in the electro-magnetic (EM) spectrum. This is the fundamental con- 
cept of electronic warfare and such importance has been placed on this mission area that 
the following maxim was developed; “The side that controls the electro-magnetic (EM) 
spectrum controls the outcome of any conflict in modern war or global politics’. Oper- 
ation Desert Storm was successful for many reasons, one of which was the allied forces 
superiority in the maintained use of systems operating at various wavelengths through- 
out the EM spectrum. Even with this superiority, there is much room for improvement. 
One example rs the Irmitation that cloud cover had on operations. Another is operating 
in regions with extensive aerosol content (such as were caused bv the intentronal burning 
of oil rigs). 

Many military devices operate in the visible to the far infrared (IR) wavelengths. 
Some examples are forward-looking infrared (FLIR), precision guided munitions and 
other electro-optical (EO) systems. For this reason, there is a great deal of interest in 
reliable estimation of EO instrumentation performance for vertical and slant-path ob- 
servations. For systems such as those above, this requires knowledge of the vertical 
variation of aerosol scattering and absorption in the wavelengths of operation. 

LOWTRAN 7 (Low resolutron transmittance model verston 7) calculates atmo- 
spheric scattering and absorption caused bv aerosols and molecules along a non- 
homogeneous path. Because there isn't one aerosol model that predicts scattering and 
absorption the best under all conditions, there are different aerosol models a user of 
LOWTRAN may choose. For maritime applications, the present model tis the Navy 
Aerosol Model (NAM). | NAM attempts to relate aerosol size distributions to 
meteorological parameters such as relative humidity, wind speed, and visibility, and then 
calculate the optical properties of the modeled aerosol. Most empirical models, such as 
NAM, do not predict the vertical distribution of aerosol required to extend extinction 
predictions to higher levels. For the maritime regime, this led to the development of the 
Naval Oceanic Vertical Aerosol Model (NOVAM ). 

The purpose of NOVAM (Gathman, 1989) is to determine the vertical distribution 
of the aerosol size spectra in the marine environment, used to obtain the associated op- 


tical and infrared properties along slant paths within the marine boundary laver. The 


means used to achieve this goal are to integrate various marine aerosol studies of a di- 
verse nature into a single model. This integration has led to the development of sub- 
models within NOVA M, one of which is selected based on the parameters that are input. 

This thesis will investigate the performance of NOVAM under conditions of weak 
cumulus convection. These conditions existed during the KEY90 experiment which took 
place from 2 July 1990 through 19 July 1990, near Marathon, Florida on the Florida 
Keys (Figure 1). 
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Figure 1. KEY90 experimental area around Marathon (shaded area) 


Il. BACKGROUND 


A. EXTINCTION PARAMETERS 
In this section, fundamental concepts and terms that form the basis of the 


LOWTRAN model, and led to the development of NAM and NOVAM are presented. 
Lambert, Beer and Bouguet established a relationship between attenuated radiation and 


the incident radiation for a homogeneous extinction as follows: 


T = L,exp( —bx) (1) 


where I is the attenuated radiation, 
I, is the incident radiation, 
b is the total extinction coefficient, 


x 1s the path length. 


To determine the total extinction coefficient, b, knowledge that extinction is caused by 


absorption and scattering by molecules and aerosol (dry or liquid) contributing linearly 


is applied as follows: 


b = D: E De E p» AB D (2) 
where 5,, — extinction due to molecular absorption, i 
b„, = extinction due to molecular (Rayleigh) scattering, 
b —^ extinction due to aerosol absorption, 
6, = extinction due to aerosol (Mie) scattering. 


In the visible region of the spectrum, 6,, and 6,, are the dominant extinction com- 
ponents. In the IR regions, 5,, and 5, are the primary components causing lowered 
transmittance (the ratio of I to /,). Figure 2 on page 5 shows the effect molecular ab- 
Sorption has on transmittance over a horizontal 1.82 km path for wavelengths up to 15 
um. This figure clearly demonstrates that atmospheric total extinction 1s wavelength 
dependant due to at least its b,, component. 

Equation 1 assumes a homogeneous total extinction throughout the atmosphere. 
Atmospheric studies have shown this assumption to generally be valid for horizontal 


paths (such as Figure 2), however, a model needed to be developed to predict the varı- 
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Figure 2. Atmospheric Transmittance (0 - 15 microns): Atmospheric 
Transmittance measured over a 1.82 km horizontal path at sea level 
(IIudson, 1969). 


ations of extinction in the vertical required to perform slant-path calculations. 
LOWTRAN is the model that resulted from this requirement. LOWTRAN was de- 
signed to calculate the total extinction as a function of altitude, and produce as output 


the radiance and/or the transmittance for a specified path through the atmosphere. 


B. THE LOW RESOLUTION TRANSMITTANCE MODEL (LOWTRAN) 

LOWTRAN’s radiance calculations account for contributions from atmospheric 
self-emission, solar and/or lunar radiance single scattered into the path, direct solar 
irradiance through a slant path to space and multiple scattered solar and/or self- 
emmission radiance into a path. LOWTRAN 7 (the latest version) incorporates separate 
molecular profiles for all major, as well as 13 minor and trace gases found in the at- 
mosphere. In addition, various aerosol, cloud and rain models are included to allow the 
user to select the models that most closely represent the region of interest. The user also 
has the option to select from six reference atmospheric profiles, or define a new atmo- 
spheric profile. A more extensive description of LOWTRAN may be found in ONTAR 
corp. (1990). (Kneizys et al., 1938) 


C. THE NAVY AEROSOL MODEL (NAMI) 
NAM represents one of the several empirically derived models contained in 


LOWTRAN and was designed to estimate the contribution of aerosol to the LO prop- 


agating characteristics in the marine environment. This model predicts an aerosol size 
distribution (dN. dr) within the atmospheric surface-layer. nominally for 10 meters, for 
specified wind speed (current and 24-hour averaged), visibihty, and relative humuditv 
(RIT) values. Three distinct types (modes) of aerosols are assumed to exist within the 
size distribution: 

|. A continental component contributed by a nearby land mass, 


2. ^ stationary component affected by the winds (controlled bv the 24-hour averaged 
wind speed), 


3. A fresh component caused by the current wind over the water. 


The predicted aerosol size distribution is obtained bv summing the contribution of each 


component as follows: 


B 
aN — Ay "A o ۳ 2 
له‎ = 2, - exp(- Int) - «0-۵ (3) 


where N represents the number density per unit radius of particles at radius ,r , 

4, 1s the amplitude parameter for mode 1, 

f. is the aerosol growth parameter for mode 1 (to be discussed in a later section), 

r, is the mode radius parameter. 
When a standard relative hunudity of 8095 is introduced into the above equation, /, is 
replaced by the value one, as will be shown in the relative humidity discussion. A plot 
of log(dN. dr) at 80% relative humidity versus log(r) is shown in Figure 3. In this figure, 
Al, A2, and AS are the amplitude: parameters as functions of air mass tvpe, 24-hour 
averaged wind speed, and current wind speed respectivelv; and r 1s the aerosol particle 
radius at the ambient relative humidity. The mode radius parameter, r,, is the radius that 
the maximum amplitude parameter value for each mode, i, occurs for this 8095 relative 
humiditv. (Gathman, 1983) 

NAM 's application 1s to make horizontal path calculations of aerosol extinction at 

the sea surface based on the following relationship between dN/dr and aerosol ex- 
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Figure 3. Idealized three lognorinal marine aerosol: (Gathman, 1983). 


where Q,, is the Mie efficiency coefficient for extinction which is a function of wave- 


ext 


length, causing f,,, to be wavelength dependent. 


ext 

To use NAM in slant path calculations, an assumed vertical size distribution must 
be developed for each mode based on the size distribution NAM produces at 10 meters. 
The present assumed size distribution profile has an exponential decrease with height. 
NOVAM is being developed to provide a capability for more accurate slant-path aer- 


osol extinction calculations in the maritime regime. 


D. THE MARITIME ATMOSPHERIC BOUNDARY LAYER (MABL) 

Vertical profiles of acrosol depend on the thermodynamic structure and processes 
within the MABL. A schematic representation of the MABL is shown in Figure 4. In 
thrs scheme, the cloud layer varics the most and, at times, may not be present. In cases 
Where the large scale flow prevents cloud development, the MABL terminates at the 
transition layer. In cases where very deep cloud development has taken place, the top 
of the MABL ts not well defined. Features such as frontogenesis and horizontal 
advection would cause the one-dimensional vicw of the MABL (as depicted in 
Figure 4) to be no longer valid. 


For a MABL as shown in Figure 4, the following principles apply: (Augstein, 1976) 
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Figure 4. Schematic of dominant processes in a multi-layered MABL: 0 = poten- 


tial temperature. (Augstem, 1976) 


Radiative cooling is compensated for mainly by the adiabatic warming of large- 
scale subsiding air and, to a lesser extent, by the sensible heat input from the ocean. 


Buoyancy driven mixing produces a well-mixed region above the shallow surface 
layer. 


Kinetic energy which is generated by the concurrent heating and moistening of air 
from below is partly converted to potential energy through entrainment of poten- 
tially warmer air from above through the top of the mixed laver. This assumes that 
buoyant parcels overshoot their level of equilibrium, thereby forcing mass of less 
density downward. 


The top of the mixed layer 1s determined by a balance of small-scale mass 
entrainment from below and the large-scale downward mass flow. 


Clouds begin to develop when convective plumes or bubbles penetrate to their 
condensation level. 


Il. THE NAVAL OCEANIC VERTICAL AEROSOL MODEL (NOVAM) 


A. BACKGROUND 

As described previously, NAM does not predict the correct vertical distributions 
which are required to extend extinction predictions to higher levels. NOVAM is designed 
to predict non-uniform and non-logarithmic extinction profiles observed in the MABL 
for wavelengths between .2 and 40 um. (Gathman, 1989) 

NOVAM is a combination of empirical and physical models which describe aerosol 
dynamical behavior. The current version of NOVAM incorporates the following mod- 
els: 


I. NAM, which defines the surface (10 meter) aerosol characteristics (Gathman, 
1983), 


2. Micro-meteorological models, which describe the vertical distribution of aerosol in 
the well-mixed MABL (Fairall and Davidson, 1986), 


3. A model that describes the vertical distribution of aerosol in the weak cumulus 
convection regime (described below) (Davidson and Fairall, 1986), 


4. A model describing the aerosol structure below stratus decks (Noonkester, 1985), 


A default boundary laver structure model of meteorological parameters based on 
surface observations (Gathman, 1989), 


6. A model that predicts aerosol growth based on relative humidity and the air mass 
type in the well-mixed and cloud layers (fitzgerald, 1978; Gerber, 1985). 
Extinction properties are calculated from the generated aerosol profiles using MIE 


scattering for which two basic characteristics exist: 


|l. The aerosol particle radius and the wavelength of incident energy are of the same 
order of magnitude. 


2. Scattering 1s generally concentrated in the forward direction (forward scatter). 


The following particle source and particle physics principles have been used in the 
development of NOVAM. Large sea-salt aerosol particles are produced from white 
water phenomena at the air-sea interface (Monahan et al., 1982). Smaller aerosol par- 
ticles are produced from sources such as gas-to-particle conversion or anthropologically 
generated aerosol. The concentration of anv particular aerosol 1s dependant on the 
source strength of aerosol production and on the mixing process as it relates to scalar 
contaminants. Regarding hygroscopic sea-salt aerosol, its size is very dependant on the 


relative humidity it 15 immersed in. Also, as the hvgroscopic aerosol picks up Water vapor 


from the atmosphere and grows in size, it changes its chemical composition and its index 


of refraction. 


B. NOVAM INPUTS 

Tables 1 and 2 show surface and radiosonde observation files (inputs) required by 
NOVAM. The surface observed parameters are used to perform the size distribution and 
aerosol extinction calculations at the surface and for sub-model selection, in conjunction 
with some information contained in Table 2. The cloud tvpe, present weather, and the 
zonal/seasonal category codes may all be found in separate tables as described by 
Gathman (1989). If any of these data items are unavailable for input, the model calcu- 
lates that parameter using default information which results in a reduced confidence in 


aerosol extinction predictions. 


Table 1. SURFACE OBSERVATION DATA FILE 





Position Meteorological Data 


Sea surface temperature (°C) 


Relative humidity ( RH) (?5) 
Current real wind speed (m/s) 
oe Averaged wind speed (24 h) (m’s) 
E. 















Air mass parameter (see text) 


8 | SS 
o feise TT 
Surface TR gm Pm O 10.6um) 

NE Present weather in standard code [0...99] 


Height of lowest cloud (m) 
Zonal/seasonal category [1...6] 





Most of the parameters listed in Table 1 are self explanatory, however, a couple 
warrant further examination. The air mass parameter is an indicator of the degree of 


continental contaminants and may be represented as follows: 


a.m.p. — ۱/4 + 1 


where Rn represents the atmospheric radon content (picocuries per cubic cm), 
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Table 2. CONTENTS OF THE RADIOSONDE DATA FILE PREAMBLE 
Row Column Description 
[1] | Number of actual radiosonde observations 
Surface potential temperature from radiosonde 
[3] | Surface mixing ratio from radiosonde 
Height of base of cloud laver (CB) 
Potential temperature just below CB 
Mixing ratio just below CB 
۱۱۱۱ Units midicator tortemperature (= Lar C) 
Potential temperature just above CB 
[3] | Mixing ratio just above CB 
Potential temperature just below CT 
[3] | Mixing ratio just below CT 
[1] | Units indicator for mixing ratio ( — lif g kg) 
Potential Temperature just above CT 
Mixing ratio just above CT 


! Minimum requirement for radtosonde data file. 





or as, 
a.m.p. = (9 x exp[ —1/4]) + 1 


where t equals the elapsed time (in days) it takes the current air mass to reach the point 


of observation from a distant land mass. 


An air mass parameter of | indicates a pure air mass (no contaminants). When the air 
mass parameter is less than or equal to 5, the aerosol is assumed to be made up of three 
lognormal components (1, 4, , and 4) as contained in the NAM description above. 
For air mass parameters greater than 5, indicating the presence of non-soluble aerosol. 


an additional lognormal component (^£) is assumed. This class of non-soluble aerosol 


is only present close to shore lines and has the same mode radius as r,, shown in 
Figure 3. Being assumed non-soluble, the mode 0 component of aerosol does not change 
size with changes in relative humidity. (Gathman. 1989) 

The remaining inputs to NOVAM are contained in a data file describing the vertical 
atmosphere, referred to as the radiosonde data file in Gathman, 1989. Table 2 shows the 
format for the first five rows that are contained in this file, necessary when choosing the 
weak convective model. Figure 5 graphically illustrates the position in the file by row 
and column designation. The remainder of the radiosonde file consists of radiosonde 
observation height, potential temperature, and mixing ratio entries. If radiosonde data 
is not available, a minimum of | surface entry must be input and a default profile is se- 
lected. 

After the input parameters are translated by NOVAM to predictions of mixing, 
source strengths, and size distributions, a prediction 1s provided of the extinction of EO 
energy as a function of altitude for the wavelength of interest. This output is based on 
the dry aerosol size distribution at a particular height and adjusted for relative humidity 


at that height. 


C. SUB-MODELS 

This section is intended to give a brief explanation of the different models being de- 
veloped within NOVAM. Figure 6 shows the major decision points and the track of 
flow of information in the current version. There are currently four sub-models in 


NOVAM: 


l. A mixed boundarv layer model (Fairall & Davidson, 1986; Davidson & Fairall, 
1986) where the boundary laver depth must be less than 3 km and not stable, 


2. A stratus model for winds less than 5 m's model (Noonkester, 1985). Here. the size 
distribution is not represented bv a series of lognorinals (as in Figure 3), there is 
no requirement for calculations to fit NAM predictions at. 10. meters, and the 
wavelengths available for NOVAM calculation are from l to 11 gm. The limita- 
tion on wavelengths 1s due to optical calculations using a simplification of Mie 
scattering parameters. 


3. A weak convection model (Davidson & Fairall, 1986) used when scattered cumulus 
clouds are present, cloud tops do not exceed 3 km, and a well-mixed laver is present 
below the cloud base, 


4. A default profile (Gathman, 1989) which is used when there 1s no vertical potential 
temperature and mixing ratio information. 


A stratus model for winds greater than 5 m/s and a deep convection model for cloud tops 
zu 


extending above 3000 m are not vet supported. 
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(Numbers in parenthesis refer to matrix location) 


Figure 5. Stylized profile definitions for use in preamble: (Gathman, 1989). 


NOVAM sub-model selection is based on the input parameters describing the ver- 
tical stratification (thermal stability, presence of an inversion, and the inversion height), 
cloud cover, cloud type, wind speed, and the requested wavelength for extinction calcu- 
lations. ۱ 

|. The default profile generator 

The default profile generator uses the surface meteorological measurements to 
develop profiles of temperature and relative humidity, generating values for every 0 
meters of altitude. The default profile differs from a radiosonde derived profile by having 
the vertical distribution of the lognormal components described by an exponential 
function of scale height. For each default altitude level, the optical calculations are based 
on the combined size distributions of the aerosol deduced from the model assumptions. 
The index of refraction and size are determined by the response (swelling) of hygroscopic 


aerosols based on the relative humidity calculated from the relative humidity profile 
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Figure 6. NOVAM flow diagram: (Gathman, 1989). 


generator (Gathman, 78). From these parameters, the aerosol extinction at each level 


can be calculated directly. 
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IV. INFLUENCES ON AEROSOL SIZE DISTRIBUTION PROFILES 
WITH WEAK CONVECTION 


A. CONCENTRATION GRADIENTS IN NOVAM 

The basis for modeling the evolution of aerosol concentration is the continuity 
equation. In the Boussinesq approximation, which. assumes a shallow planetary 
boundary laver, an incompressible fluid — 0), and negligible horizontal gradients of 
density, an aerosol component concentration, X, can be written in its simplest one- 


dimensional form as 


O) 
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where U is the mean horizontal wind, 
W is the mean Vertical wind, 
V, is the particle gravitational settling velocity, 
D is the particle molecular transport coefficient, 


WX 1s the mean vertical flux of particles due to turbulent transport. 


In this equation, X must be a conservative aerosol variable, e.g. dry size concentration, 
that is unaffected by variations in relative humidity. Terms such as IVY are unlikely to 
be known from first principles leading to applying models of the boundary layer to sim- 
plifv the process. (Davidson and Fairall, 1986) 

The Albrecht model was developed specifically for a trade-wind region (Albrecht, 
[979) and has been incorporated in NOVAM to describe the weak convective regime. 
iiteetinee-laver Structure of this regime, described by Albrecht, ts shown in Figure 7, 
In Figure 7, C represents the dry aerosol concentration with the subscripts r, b, c, and 
1 representing the surface, cloud base, the in-cloud layer, and the cloud top respectively. 


The lavers and their characteristics are summarized below, Davidson & Fairall (1986): 


Belle subeloue layer W < z < 2) where zosis the cloud base (LCD). This laver is 
characterized by linear flux profiles and “well-mixed” mean profiles. The depth of 
this laver is changed bv subsidence and entrainment. 


24. [he cumulus cloud layer (z, < z < z,), characterized by parabolic flux profiles and 
linear mean profiles. The slopes of the mean properties are determined by empirical 


models of cloud transport and entrainment properties as well as the boundary 
conditions at the top and bottom of the cloud laver. 


3. The free troposphere (non-turbulent) above the capping “trade” inversion. The 
properties of this layer are determined by advection and are usually input variables 
to the model. 


Note that the mean profiles are not continuous at interfaces z, and z,. 
Factors that influence the in-cloud dry aerosol concentration gradient, y (shown in 
Figure 7), are given by the following balance equation: 
Dy — 


—  ) ۲ ر(‎ + 


LAA, Iw AX, 
Dt E 


t(2; — 2b) (ed 


(8) 


where V.V is divergence, 
w', Is the entrainment rate of the inversion, 
t 1s the cloud relaxation time in davs, 
E is the cloud wall entrainment variable given by: 
Pe (z,— 2, )¥, t AAT ET) (9) 
AO 
where AT = SK 
۱ 
y, = the entrainment rate of @, in the laver, 


A0,, is the jump in @, at the cloud base. 


Because the layer below the cloud base is considered well-mixed, the dry aerosol con- 
centration gradient for this layer is normally considered to be 0, but a non-zero gradient 
could be predicted by the model with the specified parameters. (Davidson and Fairall, 
1986) 

Some assumptions were made to simplify calculations when incorporating the 
Albrecht model into NOVAM. The Albrecht model is dynamic, predicting the temporal 
evolution of the boundary laver. NOVAM seeks to characterize the average or typical 
vertical structure, therefore all time derivatives are neglected in the above equations. 
Additionally, the following assumptions are used to simplify calculations: (Davidson & 
Fairall, 1986) 


I. In the mixed layer, C = C 


r 
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Figure 7. Schematic diagram of aerosol concentration: | Dry aerosol concentration 
based on maritime and continental sources for a weak cumulus con- 


vection boundary layer structure. (Gathman, 1989) 


2. In the cloud layer, C = C? -- y(z —z,) 

3. C,1s given by a standard wind speed dependent model 

4. At the cloud top, assume: 
a. C? — 0 if C represents a locally generated sea-salt mode 
b. C+ ts given by climatology 

a =G eC 


These assumptions lead to the following equation for the vertical gradient from 


cloud base to cloud top: 


(Gas Cn 


en: u 


where nts a dimensionless entrainment parameter at the cloud top given bv: 


a A (11) 


In this form, the entrainment rate at the cloud base and the cloud wall entrainment are 
not important. These simplifications may reduce the accuracy in the vertical represen- 
tation of dry aerosol concentration, but this predicted vertical structure is a large first 
step in improving the exponential model currently in use. With mixing and entrainment 
controlled by the potential temperature and relative humidity profiles, a more accurate 


representation of the transport properties in the sub-cloud and cloud layers is obtained. 


B. RELATIVE HUMIDITY EFFECTS IN NOVAM 

Since aerosol particles attenuate energy in an amount proportional to their cross 
sectional area, any model requiring aerosol extinction calculations must take into ac- 
count the effect of relative humiditv on aerosol growth. Vertical changes in relative 
humidity become important when requiring slant-path aerosol extinction calculations for 
the weak convective region where relative humidity 1s expected to increase with height 
up to the cloud laver, then decrease graduallv to top of the cloud laver. Figure 8 shows 
the effect relative hunudity has on particle radius (represented by the ratio of the aerosol 
particle radius at the indicated relative humidity (r) versus the radius of the same aerosol 
particle at 095 relative hunudity (rj). Fitzgerald found that pure NaCl particles experi- 
enced a sudden increase in size at some critical value of relative humidity (between 70 
and 76% relative humidity as represented by the hatched area in Figure $) whereas na- 
turally occurring maritime aerosol particles show continuous smooth growth with 
changing relative humidity. This may be explained by the gradual deliquescence of sol- 
uble material which 1s present as a mixture of different salts. 

Based on the information shown in Figure 8, Fitzgerald developed the following 


relationship: 





S) y ae 
=| Vie ag 
where y represents the air-mass characteristic and S represents the saturation ratio (for 
the purposes of this discussion, S can be assumed to be equał to RH / 100). 
The relative humidity dependence in both NAM and NOVAM appears in the form 
of f; in the equation for size distribution (Equation 3). As can be seen in Figure $, when 


the relative humidity increases, deviation of the measured ratios from any theoretical 
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Figure 8. Growth curves for marine aerosols: Theoretical and calculated growth 


curves for NaCl particles and natural marine acrosols (Fitzgerald, 1978). 


curve which is based on dry particle size becomes greater. Because relative humidity is 
generally above 50% throughout the MABL, a reference relative humidity of 80% was 
selected to minimize the error inherent in the theoretical aerosol growth calculation, and 
this reference relative humidity was used as a basis in computing f. An equation for f 


Is generated by dividing equation 12 by the ratio of r(.8) to r, yielding (Gerber, 1985): 


fi = راز‎ EST ES, u 


where C7, C8 are conditional constants based on the aerosol mode component (recall 
Figure 3 on page 7) as displayed in Table 3. To incorporate the fourth mode (mode 0) 


into Equation 3, simply change the summation to include modes O through 3. 
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Table 3. GROWTH FACTOR CONSTANTS AND VALIDITY RANGE 


Aerosol Model Mode CT CS: Range of 
Composition Component Validitv 
(1) 


IET | — 31 — it fst [Riv 99% — 
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V. THE KEY90 EXPERIMENT 


A. GOAL 

Most observational studies within the development of NOVAM have been accom- 
plished off the California coast. Initial experiments involving the testing of NOVAM 
were performed in the same environment that the model was developed in. As expected, 
the model verifies reasonably well in that environment. Since NOVAM is intended to 
be a global maritime model, the model had to be tested in other geographical scenarios. 

The KEY90 experiment was the first attempt to test NOVAM’s performance in a 
tropic-like environment where shallow or deep convection exists. The shaded region in 
Figure 1, known as Marathon, was the region of data collection during KEY90. 
Insolation heating of the relatively shallow waters between Florida and Cuba causes 
convection to be the major meteorological process taking place near the air-sea interface 
in this region. Also, at Marathon data could be obtained bv both boat and aircraft while 


being away from land influences and major continental effects on the data. 


B. MEASUREMENTS 

Simultaneous measurements were made of parameters that the model requires as 
inputs and of extinction profiles at various wavelengths. The extinction profile meas- 
urements are then compared to the results obtained by running NOVAM at wavelengths 
corresponding to those measured. The extinction profile measurements provided the 
standard to judge NOVAM extinction profile predictions. Although this procedure 
seems fairly straight forward, it is noted that it is difficult to provide an accurate estimate 
of either the measured extinction profile or of the quality of the meteorological data used 
as inputs to NOVAM (e.g. surface IR extinction œ 10.6 um and a.m.p.). 

The data measurements collected during KEY90 were made from shore, from a boat, 
and from two aircraft. Both the atmospheric surface layer and the MABL structures 
were probed in detail. Figure 9 shows the instrumentation involved in the experiment. 
For each instrument, an asterisk marks the days of operation. Instruments are grouped 
according to the information provided (ie. air temperature, sea surface temperature, 
relative hunudity, wind speed, wind direction, atmospheric pressure, visibility, vertical 
profiles, Lidar information, and other data). Table 4 gives a summary of the time and 


location that boat and aircraft atmospheric vertical profiles were taken. For easy refer- 
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Figure 9, Summary of instrumentation used during KEY90: (Gathman, 1991) 


ence, the difference of time (in minutes) and in location (in nautical miles) 1s provided 


in the last two columns. 
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a 4. SOURCES FOR VERTICAL PROFILES 
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It should be noted that in the process of collecting data, areas of thunderstorm ac- 


tivity were avoided. This mav have caused parameters such as average wind to be lower 


than what may have actually occurred. 


C. ATMOSPHERIC SYNOPTIC SCALE 

The major synoptic flow patterns which influenced conditions in the Marathon re- 
gion from 9 to 19 July are shown in Figure 10. For the period of 9 - 10 July, surface 
winds at Marathon were generally from the east to south-east with speeds of approxi- 
matelv 3 m/s. This was a result of the influence of high pressure centered off the east 
coast of the U.S. (Figure l0a). The high pressure system continued to control the 
Marathon region winds on the 11th and 12th of July, however, short wave troughs were 
analyzed as passing through Marathon on 11 July. 

As the high pressure system moved to the north-east, winds at Marathon shifted 


from east-southeast to a more southerly direction (Figure 10b). Wind speeds were be- 
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Figure 10. Synoptic pattern in Marathon area: Surface pressure (mb) and 700 
mb heights (m * 10-! ) from NMC analysis for 1200 UTC is provided 
for the following days: a)9 July b)12 July. 
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Synoptic pattern in Marathon area (cont.): 
700 mb heights (m * 10) from NMC analysis for 1200 UTC is pro- 
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tween 3 and 5 m/s during this period. By 13 July, a low pressure system which was ap- 
proaching the eastern U.S. had an influence on winds at Marathon. Wind flow remained 
southerly, but wind speed decreased slightly due to weaker pressure gradients over the 
area. On 14 July, the low pressure center had deepened and moved to the east 
(Figure Ila). Surface wind flow over Marathon remained southerly with wind speeds 
increasing to 6 - 8 m/s. The cold front associated with the low pressure system was 
moving very slowly to the east. By late UTC on this day, it was evident that the frontal 
system was losing its upper level support, and the high pressure evident at 700 mb near 
the Florida Keys in the figure continued to build. 

From 15 to 19 July, high pressure once again dominated the flow pattern in the 
Marathon region (Figure 11b). As the high pressure center moved to the north-east, the 
pressure gradient at all levels weakened. Surface winds were from the east to south-east 
throughout this period with wind speed decreasing to approximately 5 m.s bv the I9th. 
A series of short-wave troughs were analvzed as passing through the Marathon region 
on 17 July. 

Although all observations during KEY90 fell within the realm of weak convection 
defined by NOVAM, there were examples of deep convecton in the form of 
thunderstorms. Thunderstorm activity was observed in the Marathon region on 14, 16, 
17, and 19 Julv. 


VI. NOVAM USAGE, PROCEDURE AND RESULTS 


A. PROCEDURE 
1. Surface File Generation 
Each radiosonde and flight spiral profile requires corresponding surface file. 
Figure 9 shows that there were manv choices for most surface file parameters. Because 
of consistent performance during the entire experiment, WeatherPak measurements were 
used for surface laver air temperatures, wind speeds, and relative humidities. When thrs 
information was not available on 10 and 11 July, sling psvchrometer measurements were 
used for air temperature and relative humidity while a hand held device was used for 
wind speed. Bucket thermometer measurements of sea surface temperature were used. 
NRL’s radon counter was used to establish the air mass parameter. The NOSC aircraft 
aerosol spectra was used to obtain the 10.6 um aerosol extinction at the surface. The 
altitude of this reading was usually 30 meters above the surface. [he optical visibility 
value was calculated by the model, based on the above parameters. 
2. Radiosonde Preamble Generation 
The general shapes of the potential temperature and mixing ratio (for flight) or 
specific humidity (for radiosonde) profiles were examined together when generating the 
radiosonde preamble. The first step was to establish heights of the first and second inv- 
ersions (z, and z, rn Figure 5). [n gencral, a well-mixed laver was easy to establish below 
the first inversion, where potential temperature and mixing ratio remained constant. The 
more difficult task was determining a height for the second inversion. [n most cases. 
there was no clear boundary (jump) to define the top of the cloud laver, especiallv in the 
potential temperature profiles. As a result, relative humiditv profiles and observations 
of slope changes in the potential temperature and mixing ratio profiles were considered 
subjectively in making the determination. 
3. NOVAMI Runs 
NOVAM runs for both the aircraft and radiosonde vertical profiles used the 
same surface file, so only one surface file was created for each row presented in 
Table 4. Information in the surface file is displayed above the radiosonde vertical pro- 
file panels in Figure 14 through Figure 26. 
The jumps occurring between straight line segments of the radiosonde and air- 


craft profiles are important run parameters. The predicted amount of entrainment oc- 


curring at the cloud top (second inversion) is dependent on the difference (Jumps) in the 
potential temperatures and mixing ratios at the level. The smaller the difference, the 
more entrainment is occurring between the cloud laver and the layer above the second 
inversion. The wavelength of interest for each model run was an input parameter since 
aerosol extinction varies with wavelength. All the results presented in this thesis are for 
a wavelength of 3.5 um, even though aerosol extinction profiles were obtained by NOSC 
aircraft at .53 ym, 1.06 um, 3.5 um, and 10.6 um. 

A NOVAM output file was generated for each vertical profile presented in Table 
4. A separate run was also made to generate a limiting profile. For each surface file, the 
NOVAM default vertical profile was used with the relative humidity at each level set at 
98%. When aerosol extinction measurements in the lower portion of the MABL exceed 
this value, precipitation is believed to be a possible cause, and the assumption of 
NOVAM that no precipitation exists is violated. The assumption of the default profile 
generator that each mode of the size distribution decreases exponentially with height 
will cause a difference in NOVAM versus limiting profile size distributions, increasing 
as the height of the well-mixed layer increases. 

The importance of the relative humidity to aerosol extinction can be shown by 
comparing NOVAM profiles to limiting profiles. At a wavelength of 3.5 pn, an increase 
in relative humidity from the 80% typically observed at the surface to a 98% relative 
humidity (as used in the limiting profiles) results in an approximate order of magnitude 
increase in aerosol extinction prediction (comparison of aerosol extinction at surface 
between NOVAM predictions based on vertical profiles and NOVAM limiting profiles 
in Figure 14 through Figure 26). At the times of the vertical profiles, relative humidities 
above 90% were rarely seen, therefore a comparison could not be made of how well 


NOVAM actually performed at the higher relative humidities. 


B. RESULTS 

Simple comparisons between the flight and radiosonde potential temperature, mix- 
ing ratio, and aerosol extinction profiles reveal extensive variability in space and time. 
Results from 19 Julv (Figure 26) illustrate this point the best. From Table 4, the ap- 
proximate horizontal distance between vertical profiles was five nautical miles, while the 
time difference was negligible. Analyses to determine the cloud-laver boundaries (defined 
by the first and second inversion) vielded a 150 meter difference in height of the cloud 
base and a 700 meter difference in the height that the second inversion occurred. The 


fact that there was a substantial spatial variability is further corroborated by the Lidar 
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extinction profile for 14 July shown in Figure 27, which clearly shows the spatial vari- 
ation in atmospheric structure. As discussed above, stant-path calculations of 
transmittance currently assume aerosol extinction to be horizontally uniform. This 
condition does not appear to exist in the KEY90 region. A predicted aerosol extinction 
profile generated by NOVAM based on one Vertical profile will most likely not be rep- 
resentative of the average vertical structure over the path of interest. In order for 
NOVAM to provide the required predictions, either the mean cloud base and cloud top 
levels over the path need to be established, or a technique needs to be developed to 
minimize the importance of the analvzed inversion levels. 

All NOSC aircraft aerosol extinction observations between 9 - 19 July were plotted 
simultaneously (Figure 12) to examine how well NOVAM predictions characterized 
aerosol extinction. A subjectively determined mean profile was generated from the plots 
in Figure 12, and will be referred to as the flight composite. To be considered successful, 
Separate NOVAM predictions should compare favorably to the flight composite under 
most circumstances. Due to space and tme differences between the flight profiles and 
the radiosonde profiles, onlv NOVAM results from flight profiles were used in the com- 
parisons with the observed aerosol extinction profiles. When considering only the top 


row of Figure 14 through Figure 26, the following features can be seen: 


l. NOVAM aerosol extinction predictions above the second inversion were extremely 
poor. 


2. Below the second inversion, NOVAM predictions compared favorably with the 
flight composite except in the cases where aerosol extinction dropped off noticeably 
(an order of magnitude difference between the surface and cloud base observations) 
below the cloud base (see top panel Figure I7c, Figure 20c, Figure 24c, and Fig- 
ure +25), 

The assumption that aerosol extinction above the second inversion is essentially 
zero when the air mass parameter is calculated to be approximately equal to one does 
not appear to be a good one for this region. Aerosol extinction was determined to be 
significant if it ıs within one order of magnitude of the aerosol extinction that NOVAM 
calculates at the surface for the wavelength of interest. Using this as a criteria, the 12, 
13, 14, 17, and 18 July NOVAM flight profiles all showed significant aerosol extinction 
contributions above the second inversion. 

One possible explanation for the poor performance of NOVAM above the second 
inversion is that the weak or shallow convection assumptions are not effective here ,i.e. 
this may be a case of deep convection. Thunderstorm activity and rain were frequentlv 


observed in Marathon, and multiple cloud layers were not uncommon. The hypothesis 
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Figure 12. Flight extinction profiles with flight composite: 


is that relatively small regions of deep convective activity result in an increased pro- 
duction of acrosol at the surface and transports aerosol to higher levels. As these deep 
convective cells dissipate, large amounts of aerosol are still being transported horizon- 
tally throughout the region, which would account for the aerosol present above the 
second inversion in many cases. 

Another possible explanation may be that the atmosphere above the second inver- 


sion commonly contains a significant number of small radii aerosol particles. These 
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Figure 13. | NOVAM prediction profiles with NOVAM composite: 


particles may originate either from continental sources or from the ocean surface 1n the 
form of DMS. These DMS particles would be carried to the upper atmosphere by 
convective cells and transported by large scale circulation. Particle processes associated 
with DMS origins arc not well understood. (Durkee et al., 1990) 

If the predicted aerosol extinctions above the second inversion are ignored, and a 
simultancous plot of NOVAM predictions based on flight vertical profiles is made, a 


NOVAM composite can be obtained (ligure 13) in the same manner as the flight 


>] 


composite described above. The flight and NOVAM composites turned out to be nearly 
identical (Figure 14c). There was not much variation in wind speed at Marathon and 
no frontal svstems passed through during the KEY90 experiment, therefore it 1s likely 
that aerosol generation and transport mechanisms remained fairly constant during 
KEY90. This would explain why composite mean profiles performed as well as NOVAM 
predictions below the second inversion. It is possible that mean profiles, such as the 
composites generated above, can characterize the typical aerosol extinction (at a 3.5 
um wavelength) for the Marathon region during July. This approach leads to the region 
being described empirically, therefore the predictions are specific to one region during 
one season. The required data base for describing all ocean regions for each season 
would be enormous and the time and expense involved would be substantial. 

The fact that a simple exponential curve developed from composite mean profiles 
performed much better than NOVAM above the evaluated second mversion indicates a 
need for further research on aerosol size distributions in the upper atmosphere. If the 
Structure observed in KEY90 is in fact due to deep convective processes, a rewsiem 
would have to be made to the sub-model selection process in NOVAM. Also, instead 
of the weak-convective assumption that the mode 2 and 3 components drop off to 0 
above the second inversion, an exponential function could be used to describe the aer- 
osol size distributions for these modes. ۱ 

It is believed that the observations where aerosol extinction decreased noticeably 
below the first inversion are attributed to the combination of increased aerosol pro- 
duction at a time prior to the observation, possibly caused by thunderstorm activity, and 
scavenging in the lower layer due to rain, creating a deficit of aerosol below the cloud 
layer. If future experiments reveal that rain scavenging accounts for reduced aerosol ex- 
tinction below the cloud layer, the present weather surface input file parameter is cur- 
rently the only means that this information can be relayed to NOVAM. Because the 
aerosol extinction predicted at the surface is tied to surface observations of aerosol ex- 


tinction, the task of modeling this phenomenon would still remain. 
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Figure 14. 9 July results: The top panels are based on flight vertical profiles 
observed 09 July, 2120 UTC. Bottom panels are based on radiosonde 
profiles for 09 July, 2200 UTC. Panels labeled a) represent temperature 
and mixing ratio (flight) or specific humidity (sonde) versus altitude, 
b) plots relative humidity versus altitude, and c) shows aerosol ex- 
tinction at 3.5 yen versus altitude. Refer to text for a description of the 


profiles appearing in panel c). 
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13 July results 
observed 13 July, 1330 UTC. Bottom panels 
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Figure 14. 
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same as in Figure 14. 
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same as in Figure 14. 
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same as in Figure 14. 
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Figure 23. 
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same asin Figure 14 
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on radiosonde profiles for 18 July, 1459 UTC. 
Figure ۰ 
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Figure 27. — Logarithmically amplified NRL Lidar extinction profiles: Coded in 
false color, horizontal axis is distance covered by NRL aircraft from 
1107 - 1126 UTC recorded on 14 July 1990. Vertical axis is altitude 
above sea level. Colors represent aerosol extinction in k! as given by 


the bar code on the right hand side... (Hooper, 1991) 
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VII. SUMMARY, CONCLUSION AND RECOMMENDATIONS 


An experiment was conducted in the Florida Keys from 2 - 19 July 1990 to test the 
performance of the Naval Oceanic etica] Aerosol Model (NOVAM) in a tropical, 
weak-convective regime. Meteorological data collected by aircraft and boat was used to 
generate the surface and vertical profile information files required by NOVAM. Using 
this information, NOVAM predicts the aerosol extinction (km!) for a vertical cross- 
section of the atmosphere. When NOVAM is incorporated into the LOWTRAN model, 
which also accounts for molecular extinction, slant-path calculations for transmittance 
may be performed. 

Aircraft-observed aerosol extinction profiles were obtained in order to verify 
NOVAM predictions of aerosol extinction. Comparisons between observed and 
NOVAM aerosol extinction profiles revealed major deviations above the cloud top. 
From the surface to the top of the cloud laver, NOVAM generally did an excellent job 
in predicting profile shape, with the magnitude of aerosol extinction tied to the ex- 
tinction matched at the surface. In a few cases, observed extinction increased more rap- 
idly than NOVAM predicted extinction from the surface to the base of the cloud laver. 
This is attributed to rain scavenging associated with thunderstorm activity in the arca. 
Comparison between different aerosol extinction profiles over time and space revealed 
much variation, which Lidar profiles of atmospheric structure verified. Due to these 
Variations, many vertical soundings of the atmosphere would have to be obtained in or- 
der to characterize a mean atmosphere for the region of interest. Thunderstorm activity, 
multiple cloud-layers, and the spatial variation in the atmospheric structure have led to 
a hypothesis that deep-convection was responsible for the major differences between 
observed and predicted aerosol extinction profiles. If this is the case, a simple modifica- 
tion to the weak-convective niodel may be made to apply to a deep-convection model. 

The results described above, clearly show that NOVAM is not vet ready to be used 
operationally in regions where convection is the dominant process. The following are 
recommendations for consideration in the improvement and development of NOVAM: 


1. In order to eliminate the bias in wind speed measurements caused bv avoiding areas 
of storm activitv, wind measurements from buovs placed strategicallv in the area 
of operations would be invaluable in generating the 24 hour average winds. 
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Conduct experunents in thc MARATIION area and other areas where weak 
convective conditions are expected. The data and results from these experiments 
could then be compared with the data and results from the KEY90 experiment. 


Continue to explore the possibility that some cases that are currently classified as 
weak convection may in fact be examples of deep convection. If this hypothesis is 
correct, the parameters used to differentiate between shallow and deep convection 
need to be modified. 


Explore the possibility of using satellite information to improve model predictions. 
One example would be to use the techniques established by Durkee, et al. (1986; 
1990) to determine variations in optical depth (total extinction) over a horizontal 
path which could then be incorporated into NOVAM. 


. A comparison needs to be performed between what the current model 1s producing 
(LOWTRAN,. NAM) against what NOVAM, incorporated in LOWTRAN, predicts 
(transmittance predictions versus measured would be extremely useful). 
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